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ABSTRACT 

We present a statistical analysis of the photometric properties and spatial distribution of more than 
2,800 Mgll absorbers with 0.37 < z < 1 and rest equivalent width Vl/o(A2796) > O.SA detected in 
SDSS quasar spectra. Using an improved image stacking technique, we measure the cross-correlation 
between Mgll gas and light (in the g, r, i and z-bands) from 10 to 200 kpc and infer the light-weighted 
impact parameter distribution of Mgll absorbers. Such a quantity is well described by a power- law 
with an index that strongly depends on absorption rest equivalent width Wq, ranging from ^ — 1 for 
Wo < lA to 2 for Wo > 1.5A. 

At redshift 0.37 < Zabs < 0.55, we find the average luminosity enclosed within 100 kpc around Mgll 
absorbers to be Mg = —20.65 ± 0.11 mag, which is ~ 0.5 L*. The global luminosity-weighted colors 
are typical of present-day intermediate type galaxies. We then investigate these colors as a function of 
Mgll rest equivalent width and find that they follow the track between spiral and elliptical galaxies in 
color space: while the light of weaker absorbers originates mostly from red passive galaxies, stronger 
systems display the colors of blue star-forming galaxies. Based on these observations, we argue that 
the origin of strong Mgll absorber systems might be better explained by models of metal-enriched gas 
outflows from star-forming/bursting galaxies. 

Our analysis does not show any redshift dependence for both impact parameter and rest-frame 
colors up to z = 1. However, we do observe a brightening of the absorbers related light at high 
redshift (~ 50% from Zabs ~ 0.4 to 1). We argue that Mgll absorbers are a phenomenon typical of 
a given evolutionary phase that more massive galaxies experience earlier than less massive ones, in a 
downsizing fashion. 

This analysis provides us with robust and quantitative constraints of interest for further modeling 
of the gas distribution around galaxies. As a side product we also show that the stacking technique 
allows us to detect the light of quasar hosts and their environment. 

Subject headings: quasars: absorption lines — galaxies: evolution — galaxies: halos — galaxies: 
photometry — galaxies: statistics — techniques: photometric 



1. INTRODUCTION 

Quasar absorption lines provide us with a unique tool 
to probe the gas content in the Universe with a sensitiv- 
ity that does not depend on redshift. In order to use this 
information for investigating galaxy formation and evolu- 
tion, characterizing the underlying population of galaxies 
seen in absorption is an important task to achieve. 

In 1969, Bahcall & Spitzer suggested that metal ab- 
sorption lines seen in quasar spectra are induced by large 
halos of gas around galaxies extending up to a hundred 
kpc. Since then, the study of quasar absorption lines 
has become an active field of investigation, but funda- 
mental questions regarding the nature of this gas and its 
distribution around galaxies still remain. 

Studies of absorber-galaxy connections began with 
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the investigation of Mgll systems. Such a choice re- 
sults from a practical constraint: among the dominant 
ions in HI gas, Mgll has the longest-wavelength reso- 
nance transition at the doublet AA2796.35, 2803.53. Sev- 
eral studies showed that Mgll arises in gas spanning 
mor e than five decades of neutral hydrogen column den- 
sitv (iBereeron fc Stasinskalll986l : [Steidel fc Sargentlll992t 
IChurchill et al.l [2000( 1. and at a similar time the first 
galaxies responsi ble for some me t al absorption features 
were ide ntified (iBergeronI 119861: iBergeron et al.l Il987t 
ICristianilll987t IBergeron fc Boiss ^ 119911 ). iSteidel et al.l 
(I1994D then built a sample of 58 identified Mgll sys- 
tems and found that various types of galaxies with L* 
luminosities exhibit Mgll absorption up to a radius of 
R ~ bOh^^ kpc. For a review of results regarding 
the connection between Mgll absorbers and galaxies see 
IChurchill et al.l (l2005a^. It is also interesting to mention 
that iBowen et a l. (2006) recently reported the detection 
of Mgll absorption around quasars. 

While the connection between galaxies and metal ab- 
sorbers has been clearly established, the origin of the ab- 
sorbing gas and its properties remains unclear. Whether 
these absorption lines trace gas being accreted by a 
galaxy or outflowing from it is still a matter of debate. 
Recently the SDSS database has provided us with sig- 
nificantly larger samples of absorbers (by two orders of 
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magnitude) , allowing us to accurately measure a number 
of statistical properties such as the redshift distribution, 
dN/dz, and the distribution of absorber rest equivalent 
widths, dN/dW n , which is now measured over more tha n 
a decade in Wq |Nestor et alJl2005HProchter et al.l[200l . 
Unfortunately, physical models or numerical simulations 
reproducing these observations are still lacking. While 
efforts on the theoretical side are needed, additional sta- 
tistical observational constraints are required, especially 
to characterize the spatial distribution of the gas around 
galaxies and their corresponding optical properties. 

Previous attempts to constrain such properties of Mgll 
absorbers have relied on deep imaging and spectroscopic 
follow-up of (arbitrarily faint) galaxies in the quasar field 
in order to identify a potential galaxy responsible for the 
absorption seen in the quasar spectrum. Such expen- 
sive studies have been limited to samples of a few tens 
of quasars. Moreover recent results have shown discrep- 
ancies in the inferred distributions of ab sorber impact 
parameters. Whereas ISteidel et al.l ()1994f ) claimed that 
~ L* galaxies are surr ounded by Mgll halos extending 
up to i? ~ 50/1-1 i^p^,^ IChurchill etall (|2005aD recently 
showed that a more systematic analysis reveals the same 
type of absorption up to i? ~ 100h~^ kpc in a number 
of cases. The different results apparently arise from dif- 
ferent ways of preselecting the candidate galaxies for the 
spectroscopic follow-up. Furthermore, it is interesting to 
note that the recen t detections of Mgll around quasars 
(jBowen et al.l 1200(1 ) reach similar distances. Estimating 
the distribution of absorber impact parameters without 
the use of any assumption is therefore needed. 

It should be emphasized that the identification of a 
galaxy responsible for a given absorption feature is in 
general not a well defined process. Indeed it is always 
possible that an object fainter than the limiting magni- 
tude, such as a dwarf galaxy, gives rise to the absorbing 
gas, or that more than one galaxy has a redshift consis- 
tent with that of the absorber system which is likely to 
happen in group environments. A well-defined quantity 
that can be obtained from observations and modeled the- 
oretically is the cross-correlation between absorbers and 
the measured light within a given radius: (1^0 L){9). 

In IZibetti et al.l (|2005a[ ) we proposed a new approach 
to measure the systematic photometric properties of 
large samples of absorbing systems, which combines 
statistical analysis of both spectroscopic and imaging 
datasets of the Sloan Digital Sky Survey (hereafter SDSS 
I York et aP 12000) . As Mgll lines can only be detected 
in SDSS spectra at z > 0.4, the absorbing galaxies are 
only marginally detected in individual SDSS images at 
the lowest redshifts, and below the noise level other- 
wise. However, they do produce systematic light excesses 
around the background QSOs. Such excesses can be de- 
tected with a statistical analysis consisting of stacking 
a large number of absorbed QSO images. After having 
demo nstrated the feasibili ty of this technique with EDR 
data (jZibetti et al.ll2005af l. we now extend our analysis 
to the SDSS DR4 sample, i.e., a dataset about twenty 
times larger. 

Stacking analysis has the intrinsic limitation of only 
probing the mean or median value of a quantity for which 
the underlying distribution is not necessarily known. 
However, it can lead to detections of low-level signals not 
detectable otherwise. Moreover, one can also investigate 



the dependence of an observed quantity O as a func- 
tion of a set of parameters that describe a given sample. 
Some properties of the underlying distribution giving rise 
to < O > can therefore be inferred. 

Stacking approaches have provided us with valu- 
able results in a n u mber of areas. For example, 
iBartelmann fc White! ()2003D have demonstrated that 
stacking of ROSAT All-Sky Survey X-ray images of high 
redshift clusters detected in the S PSS can be used to de - 
rive their mean X-ray properties. IZibetti et al.l ()2005bD 
detected low levels of intracluster l ight b y stacking SDSS 
cluster images, while iLTn fc Mohrl (|2004f ) measured their 
infrared (IR) luminosity by coaddin g 2MA SS data. In 
the context of galaxies, iHogg et al] ()1997D constrained 
the IR signal fr om faint galaxies using stacked Keck data. 
Similarly, Br andt et al.l (pOOl.) measured the mean X-ray 
flux of Lyman Break galaxies. IZibetti et al.l ()2004l ) char- 
acterized the very low-surface bright ness diffuse lig ht in 
galaxy halos and, more recently. Whit e et al.l ()2006f ) con- 
strained the radio properties of SDSS quasars down to 
the nanoJansky level. In spectroscopic studies, stacking 
techniques have been extensively used to look for weak 
signals. Composite spectra of SDSS quasars were re- 
cently used to detect weak absorption lines as well as dust 
reddening effects that are are well below the nois e level 
in individual spectra j Nestor et al. 2003; Mcnard 'et al.l 
I2OO5I: lYork erani200l . 

In order to emphasize the power of the stacking tech- 
nique using the SDSS, it is instructive to compare the 
sensitivity that can be reached with respect to that of 
classical pointed observations. In the present study we 
make use of SDSS imaging data which is obtained from 
a 2.5 meter telescope, using a drift-scan technique giv- 
ing rise to an exposure time of about one minute. The 
samples of quasars with and without absorbers used for 
the stacking analysis include approximately 2,500 and 
10,000 objects respectively. The careful analysis of the 
corresponding 0.5 Terabyte of imaging data allows us 
to reach equivalent exposure times of about 40 and 160 
hours respectively, i.e., 4 and 16 nights of observations 
with the SDSS telescope. Such numbers imply that, as 
long as systematic effects are small, detecting faint levels 
of light from high-redshift becomes feasible. 

In this paper we present the results of a stacking analy- 
sis aimed at characterizing the light associated with Mgll 
absorbers and investigating its properties as a function 
of absorber rest equivalent width and redshift. As we 
will show, the sensitivity achieved with this technique 
even allows us to detect the light from the high-redshift 
quasar hosts. The paper is organized as follows: in !j2|we 
present our sample of Mgll absorbers; in ^jSjwe describe 
the image processing and the technique adopted to ex- 
tract photometric information from the stack images; in 
0we analyze the surface brightness profiles of the light 
in excess around absorbed QSOs and derive the luminos- 
ity weighted impact parameter distribution of the Mgll 
systems around galaxies; |J5] is devoted to the interpre- 
tation of the integral photometry in terms of rest-frame 
SEDs and luminosity; in ^ we briefly discuss a possi- 
ble unified scenario for the interpretation of the results; 
finally, in fjTj we present a schematic summary of the pa- 
per. 

"Standard" cosmological parameters are assumed 
throughout the paper, i.e., U,tot = 1-0, f^A = 0.7, and 
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h = 0.7. 

2. THE DATA 
2.1. QSOs with Mgll absorbers 

We use the sample of Mgll abs orber system s com - 
piled with the method presented in Nestor et al. (2005') 
and b ased on SDSS DR4 data |Adelman- McCarthy et al.. 
120061 ). In this section we briefly summarize the main 
steps involved in the absorption line detection pro- 
cedure. For mor e details we refer the reader to 
iNestor et all (f2005h . AU QSO spectra from the SDSS 
DR4 database (see Richards et al. 2002 for a defini- 
tion of the SDSS QSO sample, also Schneider et al. 
2005) were analyzed, regardless of QSO magnitude. The 
continuum-normalized SDSS QSO spectra were searched 
for Mgll AA2796, 2803 doublets using an optimal extrac- 
tion method employing a Gaussian line-profile to mea- 
sure each rest equivalent width 14^0. All candidates were 
interactively checked for false detections, a satisfactory 
continuum fit, blends with absorption lines from other 
systems, and special cases. The identification of Mg II 
doublets required the detection of the A2796 line and at 
least one additional line, the most convenient being the 
A2803 doublet partner. A 5a significance level was re- 
quired for all A2796 lines, as well as a 3a significance level 
for the corresponding A2803 lines. Only systems 0.1c 
blueward of the quasar redshift and redward of Lyman 
a emission were accepted. For simplicity, systems with 
separations of less than 500 km/s were considered as a 
single system. We note that Mgll absorption lines are in 
general saturated, and in these cases no column density 
information can be directly extracted from Wo('^2796). 
Given the typical S/N of SDSS quasar spectra, most of 
the detected Mgll absorption fines have Wq > 0.8 A. 
Weaker systems can occasionally be detected in the spec- 
tra of bright quasars but are significantly less numerous. 
As our technique requires a large number of absorbers 
we will only consider the stronger ones, i.e., Wq > 0.8 A. 

In this sample, multiple absorbers detected in the same 
quasar spectrum are found in ~ 15% of the systems. In 
order to ease the interpretation of our analysis, we select 
quasars with only one strong absorber detected in their 
spectrum. Finally, we restrict our analysis to Mgll sys- 
tems with z < 1. After excluding a few tens of QSOs for 
which poor imaging data is available (mainly because of 
very bright saturated stars close by), the sample of ab- 
sorbed QSOs includes 2844 objects. Their properties are 
summarized in Figure[T] the left panel shows the redshift 
distribution of the selected sample of Mgll absorbers 
(shaded histogram) and their background quasars (solid 
line). The right panel presents the distribution of ab- 
sorber rest equivalent widths. 

2.2. Reference QSOs 

In order to use a control sample in our analysis we need 
to select a population of quasars without absorbers. A 
number of criteria must be applied to ensure the selec- 
tion of an unbiased sample: the two quasar populations 
must have the same redshift and apparent magnitude 
distributions, and must present the same absorption-line 
detectability. The selection of the reference quasar popu- 
lation is therefore defined as follows: for each quasar with 
a detected absorber with rest equivalent width Wq and 



redshift Zabs , we randomly look for a quasar without any 
detected absorber above Wo(A2796) = 0.8 A such that: 

(i) they have similar redshifts: Az < 0.1; 

(ii) their magnitudes do not differ significantly: Arrii < 
0.5; 

(Hi) the absorption line is detectable in the reference 
spectrum, i.e., the S/N of the reference spectrum 
at the corresponding wavelength is high enough to 
detect the Mgll A2796 line at the required level 
of significance. Such a requirement translates into 

WminiQSOref , Zabs) < Wq, whcrC Wmzn is the 

minimum equivalent width that can be detected 
by the line finder at the corresponding wavelength. 

The quasars selected in this way are called reference 
quasars in the following. Our selection procedure en- 
sures that if similar absorbers were present in front of 
the selected reference quasars, we would have been able 
to detect them. In order to improve the noise proper- 
ties of our analysis we select four independent reference 
quasars for each quasar with an absorber, resulting in a 
total of 11,376 objects. 

2.3. Definition of absorbers sub-samples 

The complete sample of Mgll absorbers used in this 
study spans a large range in Zats (from 0.37 to 1.00), 
where substantial wavelength shift occurs. In fact, at 
the lowest redshift, Zabs = 0.37, the 4000A break feature 
fies in the observed r-band, whereas at Zats = 1-00 it 
occurs in the observed z-band. Moreover, the distance 
modulus between Zabs — 0.37 and Zats = 1-00 increases 
by 2.6 mag, thus making low-Zahs absorbing galaxies the 
dominant sources of signal in the overall stack. There- 
fore, in the following analysis we will consider three 
Zabs bins: 0.37 < Zabs < 0.55 ( "low-redshift" sub- 
sample), 0.55 < Zabs < 0.76 ("intermediate-redshift"), 
and 0.76 < Zabs < 1-00 ( "high-redshift" ) . The bound- 
aries of this binning in Zabs are indicated by the vertical 
dashed lines in the left panel of Figure [1] Possible depen- 
dence of the photometric and spatial properties on the 
absorption strength will be investigated by further split- 
ting each Zabs subsample into three different Wo(A2796) 
subsamples: 0.8A < Wq < 1.12A ("Iow-Wq" systems), 
I.I2A < Wo < 1.58A ("intermediate- Wo" systems), and 
Wo > 1.58A ("high- Wo" systems). The boundaries of 
these bins are also shown as vertical dashed lines in the 
right panel of fig. [TJ Finally, it is important to note 
that very strong systems (Wo > 3A ) are not expected 
to contribute significantly to our results as they are sub- 
stantially less numerous than weaker ones. 

3. IMAGE PROCESSING AND STACK PHOTOMETRY 

As indicated above (see also lZibetti et al.ll2005aD . our 

approach makes use of the fact that the galaxies linked 
to an absorbing system statistically produce an excess 
of surface brightness (SB) around the absorbed QSOs 
with respect to unabsorbed ones. Such a SB excess can 
be measured to obtain both global photometric quanti- 
ties for the absorbing galaxies in different bands and the 
spatial distribution of the light of the absorbing galaxies, 
from which the impact parameter distribution of absorb- 
ing gas clouds can be derived. 
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Fig. 1. — Left panel: redshift distributions of Mgll absorbers with Wq > O.SA used in this analysis (shaded histogram), of their 
background quasars (soUd Hne) and of the complete population of reference quasars (dashed line). The latter is rescaled by a factor 4 to 
allow a direct comparison with the population of absorbed quasars. The vertical dashed lines mark the boundaries between the three bins 
in absorbers' redshifts. Right panel: distribution of absorber rest equivalent widths. The vertical dashed lines separate the three bins in 
absorption strength. 

In this section we describe the techniques that allow 
us to optimally integrate and measure the flux distribu- 
tion of all and only the galaxies that cross-correlate with 
the presence of an absorber. By "optimally" we mean 
that any source of noise is minimized. In this measure- 
ment we are faced with three main sources of noise: (i) 
the intrinsic photon noise, which is fixed by the number 
of stacked images and their depth, (ii) the "background" 
signal produced by stars and galaxies which are not corre- 
lated with the absorbers, and (Hi) the light from the QSO 
itself. While the intrinsic photon noise appears to be 
sufficiently low for a few hundred SDSS images and ab- 
sorbers up to 2 = 1 or more, the signal from background 
sources and the brightness mismatch which is allowed 
between absorbed and reference QSOs, although small, 
produce a noise which is orders of magnitude larger than 
the signal to be detected. These two sources of noise 
must be drastically reduced by applying accurate mask- 
ing and QSO/Point Spread Function (PSF) subtraction 
algorithms on each image as described in the next two 
sections i) iJ3. II and 13.21 

The stacking and subsequent photometric analysis is 
conducted simultaneously in the four bands for which 
most of the flux is expected, i.e. g, r, i, and z. We use the 
so-called SDSS "corrected frame" images (fpC's), which 
are bias subtracted and flat-field corrected, but preserve 
the original background signal. The error computation 
takes into account not only formal photometric uncer- 
tainties, but also sample variance and error correlation 
between different quantities in different bands. This is 
described in detail in 



3.1. Masking algorithm 

The goal of the masking algorithm is to mask out any 
region of the image frame where flux is contributed by 
sources that have negligible probability of being a galaxy 
at the redshift of the absorber, Zabs- As a preliminary 
step, the light of the QSO is removed from a working copy 
of the image. This is done by subtracting the so-called 



postage-stamp image of the QSO output by the SDSS 
PHDTD pipehnc (jLupton et al.ll2001f ). in which the QSO is 
optimally de-blended from all surrou nding sources. SEx- 
TRACTOR (jBertin fc Arnoutsl[l996l V.2.3.2) is then nm 
on the QSO-subtracted image to obtain a segmentation 
image and a catalogue of all sources detected la above 
the local background over a minimum area of 10 pixels. 
The SExtractor's segmentation image is used as base 
for the final mask. Only segments (i.e., groups of pixels) 
attributed to sources that are too bright to be a galaxy 
at Zabs are kept in the mask. 

For the mask thresholds, we adopt the fluxes derived 
for an unobscured, metal poor stellar population pro- 
duced in a 10 Myr long burst, observed right after the 
end of the burst itself, with rest-frame g-hand absolute 
magnitude Mp tvirsshniri = ~22.4 mag (as computed by 
the code of iBruzual fc Char loti 12003'). observed at Zats- 
The flux- limited masks are computed separately in g-, 
r- and i-band, that roughly correspond to the rest-frame 
g- and u-bands in the redshift range spanned by the ab- 
sorbers in our sample. The choice of this very "blue", 
UV-bright SED is motivated by the need to avoid the 
exclusion of UV-bright galaxies, while keeping the flux 
thresholds relatively low in order to minimize the con- 
tamination by bright interlopers. Compared to the field 
luminosity function (LF) in a similar redshift range (cf. 
'Gaba sch et"alll2004l ). our flux thresholds are ~ 1 mag 
brighter than M* in the rest-frame g-band, and ~ 1.8 
mag brighter in u. We make sure that the exclusion of 
galaxies at the very bright end of the LF does not af- 
fect our results by performing the stacking using a more 
conservative value of Afg^thrcshoid = —23.9 mag (i.e. 1.5 
mag brighter than previously assumed). While the noise 
from bright interlopers increases significantly, no system- 
atic effect results from this change in the masks. There- 
fore, all the analysis of this paper is conducted on images 
masked with Mg^throshoid = —22.4 mag. For illustration, 
we plot the apparent r-band threshold magnitude as a 
function of the absorber's redshift (solid line) in Figure 
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Fig. 2. — The threshold magnitude for the masks in r-band as 
a function of the absorber's redshift Zabs (soUd Hne), compared to 
the apparent magnitude of a "normal" (Sbc SED template) Mg = 
M* galaxy (dotted line). The horizontal long dashed line is the 
magnitude where completeness of the SDSS photometric catalog 
drops below 95% (point sources). 

[21 The dotted line displays the magnitude of a Mg = M* 
galaxy for an intermediate (Sbc) SED template. The 
choice of a UV-bright template yields flux thresholds that 
become increasingly brighter than "normal" M* galaxies 
at increasingly higher redshift. 

It is worth noting that at low redshift some absorbing 
galaxies can possibly be detected in the SDSS images: in 
fact, the SDSS photometric sample has a completeness 
larger than 95% (for point sources) at r < 22.2 mag and 
thus would allow to detect objects as faint as Mg = M* 
up to z = 0.65, and 10 times fainter than that at z = 0.37 
(an analysis of galaxy number counts will be presented in 
a forthcoming paper). However, at Zabs ^ 0.7 the detec- 
tion of individual galaxies is completely precluded and 
only the stacking technique is able to reach the required 
sensitivity. In the second step, we build a mask where 
all stars brighter than 21.0 mag in r are covered with a 
circular patch of radius R = 1.2 x 025 + 3 (in pixels), 
where 025 is the median of the isophotal semi-major axis 
at 25 mag arcsec"^ in g, r, and i. For saturated stars, 
is replaced by three times its value, in order to avoid 
the extended halos of scattered light that surround these 
stars. For this step we entirely rely on the SDSS-DR4 
photometric catalogs (PhotoPrimary table) in order to 
ensur e optimal star/galaxy separation (see llvezic et al.l 
[200I . 

Finally, the flux-limited masks in the three bands and 
the star mask are combined into a final mask, where a 
pixel is masked if and only if it is masked in any of the 
input masks (logical or). 

3.2. QSO's PSF subtraction 

In very high S/N images, like the stacked images, the 
spread light of a QSO can be detected out to angular 
distances of a few tens of arcsecond and still contribute 
a non-negligible fraction of the azimuthally averaged SB 
linked to the absorbers. In order to suppress this con- 
tribution with the best possible accuracy, the SB distri- 
bution of the QSO is estimated from the high S/N PSF 



of bright, unsaturated stars. More specifically, for each 
QSO we look for a star that fulfills the following require- 
ments, (i) The star must not be saturated, nor fainter 
than 17 mag (r band); stars with interpolated pixels, 
cosmic rays, or flagged for any photometry defect are re- 
jected, as well as those stars whose centers happen to 
lie less than 5 pixels apart from any mask in the frame 
(see below how PSF star frames are masked), (ii) The 
star must have been observed during the same run as 
the QSO, in the same SDSS camera column, and within 
10 arcmin angular separation from the QSO. These two 
criteria, (i) and (ii), ensure that the PSF is photometri- 
cally accurate and it is taken in observing conditions that 
are consistent with those of the QSO. Among the stars 
selected in this way, the one that best matches (Hi) the 
size (as determined by the second moment of the light 
distribution) and colors [g — r and r — i) of the QSO is 
chosen as representative of the QSO's PSF. This last se- 
lection aims at minimizing the effects of the space/time^ 
and color dependence of the PSF. In fact, it must be 
noted that QSOs are on average bluer than stars and 
therefore, have broader PSFs at any given position, on 
average. The selection at point (Hi) minimizes this effect 
and, at the same time, compensates for it by exploiting 
the spatial variation of the PSF. The stars selected to 
model the PSF have median r ~ 15.7 mag, and span the 
range between 13.4 and 17 mag. The median difference 
in color between QSOs and the corresponding stars is 
—0.28 mag (rms 0.21) in 5 — r and —0.06 mag (rms 0.16) 
in r — i. The distribution of difference in size is centered 
at 0, with 0.16 pixels rms. 

The masking algorithm described above ( ijS.ip is then 
applied to each selected reference star; however, no lower 
flux threshold is applied and all SExtracted sources are 
retained in the masks. The subtraction of the PSF frame 
from the QSO frame is only possible for pixels that are 
unmasked in both frames. Thus, in order to limit the loss 
of effective detection area around the QSOs, we "restore" 
as many as possible masked pixels in the PSF frame in 
the following way. The PSF frame is virtually divided 
into four quadrants centered on the reference star. Each 
masked pixel is then replaced by the average of all (if 
any) unmasked pixels at axisymmetric positions in the 
other three quadrants. If none of the symmetric pixels is 
unmasked, the pixel is kept masked. This "symmetriza- 
tion" process is very effective and results in no more than 
0.05% of the area masked in 95% of the PSF frames. 

The normalization of the PSF to the QSO, necessary 
before the subtraction can be done, is computed using 
the PSF magnitudes from the SDSS-DR4. 

3.3. The stacking and photometry procedures 

Before all images are actually stacked, we rescale them 
to the same physical scale at the redshift of the absorber, 
Zabs- This is determined such that at the lowest redshift, 
Zabs = 0.37, the original SDSS pixel scale is preserved 
(0.396 arcsec pixel"^); the resulting common pixel scale 
is thus 2.016 kpc pixel"^. Pixel intensity is accordingly 
rescaled in order to conserve the total flux in any ob- 
ject. The images of reference QSOs are rescaled accord- 
ing to the redshift of the absorber in the corresponding 

As the SDSS images are acquired in drift scan mode, a position 
offset in the scan direction implies a time delay as well. 
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absorbed QSO. PSF stars are rescaled with the same ge- 
ometry as their corresponding QSO. 

As individual images contribute to the final stacked 
image only with their unmasked regions, it is mandatory 
that all images are accurately background subtracted to 
avoid spatial background fluctuation in the flnal stacked 
image that might strongly reduce the photometric accu- 
racy. A constant background level in each image (QSOs 
and stars) is computed as the median of the unmasked 
pixel values within an annulus between 400 and 500 kpc. 
Two (j-clipping iterations with ±3(T cut are performed to 
exclude the faint sources that are left unmasked by the 
masking algorithm (see above, §3.1|) . 

All images are also intensity rescaled in order to match 
a common photometric calibration (accor ding to the 
SPSS photometric standard defined by Fukugita et al.1 
119961 : ISmith et al.l I2OO20 . that includes the correction 
for Galactic extinction. More specifically, the calibrated 
pixel intensity Icaiib is derived from the raw pixel inten- 
sity in the fpC Iraw according to the following equation: 

^calih — ^raw ' « ' ru yi j 

where /20 and /2o,rcf are the pixel counts correspond- 
ing to 20 mag (observed) in the fpC original frame and 
in the final reference photometric system, respectively. 
A\ (in mag) is the galactic a ttenuation in the wo rking 
pass-band, computed from the lSchlegel et all ()1998f ) dust 
distribution map. 

The final stacked image in each band, for any given 
sample of QSOs (either absorbed or reference), is com- 
puted such that the resulting intensity at any pixel is the 
simple average of the net PSF-subtracted counts over all 
corresponding pixels which are unmasked in the individ- 
ual images: 

r _ Eili mask{i) x {Iqso,i ~ IpsF,i) 

^stacked— ly ) 

Z^j=i mask[i) 

where Istacked is the intensity of a given stacked 
pixel, lQso,i and IpsF,i are the calibrated, background- 
subtracted intensities of the QSO and PSF frames at the 
corresponding pixel in the i^^ image, and mask{i) equals 
1 if the pixel is unmasked in both the QSO and the PSF 
image, otherwise. It is worth mentioning that, at a 
given position on a stack image, only 6 to 11% of the in- 
dividual images are masked in the \ow-Zabs bin and even 
less at higher redshift. The final intensity, Istacked, is 
therefore a reliable estimate of the average value. A care- 
ful re-determination of the (residual) background level is 
done on the final stacked image in an annulus between 
350 and 500 kpc. Circular aperture photometry is then 
performed, by simply averaging the pixel intensities in a 
series of circular annuli. 

For the purpose of determining uncertainties that take 
into account not only formal photometric errors, but also 
the effects of sample variance and the correlation between 
different apertures and bands, the full covariance matrix 
is computed by means of a jackknife algorithm. For any 
given subsample of N QSOs, the stacking and photo- 
metric measurements in the four bands are repeated N 
times, leaving out a different QSO each time. The covari- 
ance matrix of the aperture photometry fluxes obtained 
from these N realizations is then multiplied by — 1 to 



give the sample covariance matrix. All the errors rela- 
tive to the photometric quantities derived from the pri- 
mary aperture fluxes, such as integrated fluxes, colors, 
and proflle shape parameterizations, are computed using 
this covariance matrix. 

Absorbed and reference QSOs are stacked and ana- 
lyzed separately as far as primary aperture photometry 
is concerned. However, all the results that will be pre- 
sented and discussed in the following sections are ob- 
tained from net (i.e., absorbed— reference QSOs) aper- 
ture photometry. For the error estimates of the derived 
quantities, absorbed and reference QSOs are treated as 
independent datasets, thus the combined covariance ma- 
trices have null cross terms. We note that this assump- 
tion is correct only as a flrst order approximation; as 
we will discuss below in t j4.1l and Appendix [Xl there are 
systematic effects in the stacked reference QSOs that de- 
pend on the redshift and magnitude of the QSO, thus 
making the absorbed and reference samples correlated. 
Neglecting such second-order error correlation results in 
slightly over-estimating errors on the flnal net photom- 
etry, especially within the seeing radius (i.e., the radius 
that corresponds to the FWHM of the PSF image). 

All the stacking and photometric measurements de- 
scribed in this subsection and in the following sections 
are performed using a dedicated package of programs and 
utilities, which are developed in C lang uage by S.Z. The 
package relies on the cf itsio libraries (|Pencelll999r ) for 
input/output interface an d makes use of a n umber of Nu- 
merical Recipes routines (Press et al.lll99^ . The analy- 
sis presented in this paper results from the processing of 
0.5 Terabyte of imaging data. 

4. THE SPATIAL DISTRIBUTION OF MGII ABSORBING 
GALAXIES 

In this Section we present an analysis of the surface 
brightness distribution of the excess light around ab- 
sorbed QSOs and infer from that the spatial distribution 
of absorbing galaxies with respect to the Mgll systems 
probed by the QSOs sight-lines. In particular, in []4T]we 
present the SB profiles extracted from the stack images, 
and parameterize their shape as a function of the absorp- 
tion rest equivalent width, Wo(A2796), and absorption 
redshift, Zabs- In ij4.2l we will show that the SB profile 
is directly related to the impact parameter distribution 
weighted by the luminosity of the absorbing galaxy. Such 
a distribution will be explicitly derived and discussed in 

El 

4.1. Stack images and surface brightness profiles 

In Figure [3] we report an example of typical final stack 
images for the absorbed QSOs (left panel) and for the 
corresponding control sample of reference QSOs (right 
panel). More specifically, this is the r-band stack for the 
low Zabs sample that includes all Wq. The two panels 
reproduce a region of 600x600 kpc^ projected at Zabs, 
centered on the QSOs. The grey level intensity maps 
the surface brightness linearly. The corresponding lev- 
els, expressed in physical units of apparent magnitudes 
per projected kpc^, are indicated by the vertical color bar 
between the two panels. The infinity level (00) indicates 
the background signal. In both panels, we note that the 
intensity of the very central pixels is zero, within the typ- 
ical background noise uncertainties. This is obtained by 
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Fig. 3. — The r-band stack images of absorbed {left panel) and corresponding reference QSOs {right panel) in the redshift range 
0.37 < Zabs < 0.55. The sample of reference quasars is four times larger which results in a lower noise level. Both images are PSF-subtracted 
and are visualized with the same color scale, as indicated by the color bar marked in units of magnitudes per kpc'^ (oo corresponds to the 
background level). The excess of light around the absorbed QSOs is apparent. 



construction, by subtracting the PSF from each individ- 
ual stacked image. The fact that these central pixels are 
zeroed to such high accuracy is indicative of the good 
quality of our PSF subtraction algorithm. What is most 
important to notice in Figure [3] is the SB excess around 
the QSOs, beyond the very central pixels. Although an 
excess is seen in both panels, it is apparent that the SB 
excess around absorbed QSOs is much more intense and 
extended than around the reference QSOs; it can be seen, 
even by eye, that the excess around the absorbed QSOs 
extends well beyond 50 kpc, whereas nothing is appar- 
ent around the reference QSOs beyond 10-15 kpc. This 
effect is even more significant if one considers that the 
reference QSO image results from 4 times as many QSOs 
as in the absorbed QSO image (see §2.2p . and therefore 
has twice as low background noise. 

The SB excess observed around the reference QSOs is 
a systematic effect which is mainly caused by the emis- 
sion of the QSO's host galaxy /environment. We ana- 
lyze the origin of this effect in greater detail in appendix 
[XI Regarding the measurement of the light of absorb- 
ing galaxies, we can treat the SB excess found around 
reference QSOs as a pure systematic additive contribu- 
tion, which is inherent to any subsample of QSOs with 
a given redshift and magnitude distribution. Therefore, 
for each subsample of absorbed QSOs, the final stack im- 
age (SB profile) of the corresponding reference QSOs is 
subtracted from the absorbers ' one to obtain a net image 
(SB profile). 

Next we analyze the SB distribution of the galax- 
ies responsible for the Mgll absorption by extracting 
one-dimensional SB profiles from circular annuli, within 
which the average SB is computed. The spacing between 
subsequent apertures is chosen to ensure adequate S /N in 
all apertures; this results in a progression that increases 
slightly faster than a geometric one. The SB profiles 
of the absorbers in all our subsamples are presented in 
Figures [4] to [21 for the g, r, i, and z band, respectively. 
Each of these figures is organized in twelve panels: three 
rows corresponding to the three Zabs-bins (from the low- 
est in the bottom row to the highest in the top one) , and 
four columns for the low-, intermediate- and high-Wo 



systems, plus the total of the Zabs-bin ("all"), from left 
to the right. The impact parameter R is reported on the 
x-axis in logarithmic scale, while the surface brightness is 
given on the y-axis in units of apparent magnitudes per 
projected kpc^. The net SB profile is shown by the black 
filled dots with error bars. Vertical error bars represent 
the jackknife estimate of the SB uncertainty computed 
as above f H3.3p . while the horizontal error bars just mark 
the width of the radial bin. The orange open circles with 
upper and lower ticks indicate the uncorrected SB deter- 
mination and relative jackknife uncertainty, prior to the 
subtraction of the systematic excess measured from the 
reference QSOs sample. It must be noted that beyond 50 
kpc no subtraction is performed (hence no orange points) 
as the reference stack images are completely consistent 
with zero at this impact parameter, and the subtraction 
would just increase the noise in the net profile. 

As a preliminary sanity check, we note that the sig- 
nal decreases from the lowest to the highest Zabs bin by 
roughly 2 mag, which is consistent with the cosmological 
dimming expected for the covered redshift range. The 
first important result that we can extract from Figures 
[31 to [7| is that there is light that cross-correlates with the 
Mgll absorbing clouds out to 100-200 kpc projected dis- 
tance. This clearly indicates that Mgll absorption sys- 
tems probe the intergalactic medium out to quite large 
galactocentric distances, and may possibly arise in com- 
plex galaxy environments. 

It must be noted that the PSF subtraction artificially 
introduces the constraint SB = within the first reso- 
lution element (corresponding to the seeing). In Figure 
[31 the two vertical short-dashed lines in grey reported in 
each panel show, at the lowest and highest Zats in tbe 
bin, the value of 1.4" which corresponds to the median 
seeing (FWHM) in SDSS images. Due to the above nor- 
malization we expect the signal to be strongly decreased 
within this radius (~ lOkpc), and significant deviations 
from the "true" SB profile to occur out to twice this dis- 
tance from the QSO 20kpc) as a result of the breadth 
of the seeing distribution in our images. 

The SB profiles between i? 20 kpc and i? = 100 
kpc in all bands are generally well approximated by a 
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Fig. 4. — The g-band surface brightness profiles of the Hght around absorbed QSOs, in different bins of redshift z^bs (indicated by the 
labels on the right-hand side) and rest-frame equivalent width of the 2796A line (labels along the top of the figure). Black filled dots with 
error bars represent the net SB, i.e., after subtracting the residual systematic contribution measured from the reference QSOs inside 50 
kpc. The orange open circles with error bars show the SB prior to the subtraction of the SB of the reference sample. The dashed green line 
is the same in all four panels of each row, and displays the best power-law fit for the sample including all Wo(A2796) in that z^^^s bin. The 
power-law fits for the three Wo(A2796) subsamples in each redshift interval are shown by the red solid lines. In each panel, vertical dotted 
lines mark the distance corresponding to the 1.4" average seeing (FWHM) at the minimum and maximum redshift in the subsample. 



single power-law with index ranging from —1 to —2.5. 
Between 10 and 20 kpc (roughly corresponding to one to 
two times the median seeing value at all redshifts) we ob- 
serve a break and a deficit with respect to the powerlaw, 
which can be mostly ascribed to the central normaliza- 
tion. However, real departures from a powerlaw profile 
in the inner 20 kpc cannot be excluded based on the 
present analysis. Also, it is worth noting that previous 
studies show that very few galaxies, if any, are associ- 
ated with Mgll absorbers at impact paramete rs < 10 
kpc (Stei del"eral1[l99l IChurchill et all l2005af) . Thus 
the signal that we miss from the central part of the pro- 
file is probably negligible. Beyond 100 kpc the photo- 
metric uncertainties are too large to derive any reliable 
information about the profile shape, although significant 
detections are obtained in several bands and bins. 



Power-law fits of the form 

S{R) = S(R) 



( = 



(3) 



are performed in the impact parameter range 20 kpc < 
R < 100 kpc by minimizing the computed using the 
full data covariance matrix. The slopes, a, for each 
absorber bin are reported in Table [T] (columns 2-5, for 
each band respectively), along with their upper and lower 
confidence limits (l-c), which are obtained from the 2- 
dimensional distribution in the S{R) — a parameter 
space. Column 6 reports the weighted average of the 
slopes in the four bands, with relative error''. The best 
fitting power-laws are also shown in the panels of Figure 

^ We note that the error that is quoted in this case is slightly 
under-estimated because we neglect that the slope determinations 
in the four bands are not completely uncorrelated. 



The light of Mgll absorbing galaxies 



9 



r-band 



0.80A^Wo<1.12A 1.12A^Wo<1.58A Wg^LSSA 



ALL Wo^O.BOA 

1 1 1 1 — I I I 1 1 iij r 



O 

Oh 




- o 



-v2 

03 



- 1^ 



A 



O 
O 



o 

cn 

1^ 



- A 



03 



CO 
-n2 

1^ 

N 

Ql 

cr 
u 

A 
O 

Qi 

CP 



100 10 



Fig. 5. — Same as Figure l4l but for the r-band. 



100 10 

R [kpc] 



10 100 



13] to [71 They are shown in red for the three EW bins^ 
and with a dashed green line for the "All"-EW pan- 
els. The fit to the "All"-EW sample is also reported 
as a dashed green line in all other panels to facilitate a 
prompt comparison between different subsamples. 

Looking at the "All"-EW samples we observe that, 
within the quoted uncertainties, all profiles in the four 
observed bands and in the three redshift bins are consis- 
tent with a common slope, a, between —1.4 and —1.6. 
There is no evidence for any redshift evolution of the 
power-law slope for the average SB distribution of ab- 
sorbers related galaxies between Zabs — 0.37 and 1. This 
is confirmed and better constrained by the weighted av- 
erage of the slopes in the four bands reported in column 6 
of Table [TJ their values are consistent within 1-a (< 0.17) 
among the three Zabs bins. A second interesting result 
that emerges from Figure |4] to [7] can be derived by com- 
paring different EW bins, especially in the most sensitive 
r and i bands for the low-Zabs systems. It is apparent that 

No red line is reported in three panels for the z band (Figure 
[TJi as the very low S/N of the profiles resulted in l-cr uncertainty 
range on the power-law slope larger than 3. 



the high-VKo Mgll systems have systematically steeper 
power-law profiles than the weaker systems. Despite the 
relatively large uncertainty of our slope measurements, 
this trend is clear in the low- Zafcs bin from the values re- 
ported in Table [H We also plot these values in the top 
panels of Figure [51 where different symbols represent the 
measurements in the four bands as a function of the ab- 
sorption strength, binned as defined in tJ2.3l The three 
columns of panels are for the three Zats bins. The large 
filled diamonds in the lower row represent the weighted 
average of the four bands, while the values for the dif- 
ferent bands are shown in the upper row. At the lowest 
Zabs the trend of steepening of the profile with absorption 
strength is unambiguously detected in r and i and statis- 
tically confirmed by the weighted averages. The power- 
law slope steepens from —1.1 in the weakest systems, to 
—2.3 in the strongest. Similar trends (in particular for 
the weighted averages) are seen at higher redshift too. 

We further investigate the systematic variation of the 
surface brightness profiles as a function of the absorption 
strength in a model independent, non-parametric way, by 
means of the first moment of the SB distribution in R. 
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Fig. 6. — Same as Figure l4l but for the i-band. 



This is defined as: 



Ri 



100 kpc 

/ S{R) dR 

10 kpc 
100 kpc 

/ S{R) R dR 

10 kpc 



Ri is thus given in kpc and can take any value between 
10 and 100 (see above in this section why 10 kpc is chosen 
as lower limit) . As we will demonstrate in the next ij4.21 
i?i can be interpreted as the luminosity- weighted average 
impact parameter of Mgll absorbing galaxies. Small val- 
ues of Ri characterize centrally concentrated light pro- 
files, while less concentrated profiles have larger The 
values of i?i computed for the four bands and their sta- 
tistical errors are reported in columns 7 to 10 of table [U 
for all bins; the weighted average among the four bands^ 
is reported in column 11. We note that i?i is remark- 
ably constant among different redshift. For "All" EW 

^ No correlation between the four bands is taken into account 
in this case, possibly resulting in underestimating the associated 
error by a few per cents. 



100 10 100 10 100 

R [kpc] 



< i?i > is 48 kpc. However, for strong systems this value 
decreases to 43 kpc (at \ow-Zabs, similar at higher red- 
shift) while it approaches 60 kpc for the Iow-Wq systems. 
As for the power-law index, we also plot these values in 
fig. [8] (bottom panels). Again, these plots confirm that 
stronger systems have more centrally concentrated pro- 
files. The effect appears statistically significant at low 
and intermediate Zabs, especially when the weighted av- 
erage of the four bands is considered. For the highest 
redshift bin the trend is much more uncertain because 
of the lower S/N, but is still consistent with the lower 
redshift samples. 



(4) 



4.2. From surface brightness profiles to the impact 
parameter distribution of Mgll systems 

In this section we provide a derivation of the link be- 
tween the surface brightness profiles obtained from the 
stacked images and the impact parameter distribution of 
the Mgll absorbing clouds. First, we demonstrate that 
if each Mgll absorber is linked to only one galaxy, then 
the SB at any radius R is proportional to the luminosity 
weighted impact parameter distribution. 
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Fig. 7. — Same as Figure |4] but for the z-band. Note that for three out of the nine Wo subsamples no power-law fits (red lines) are 
reported, because the large uncertainties make the fit meaningless. 



Let us consider initially the case of galaxies with fixed 
luminosity L, surrounded by Mgll absorbing clouds re- 
sulting in a probability of producing an absorption along 
the sight-line 



ness: 



SB{R, L) 



dNabs{R,L)-L 



(8) 



dP{R,L) = p{R,L) dA, 



(5) 



where p{R, L) is the differential probability per unit area 
at given luminosity X, which we assume to be indepen- 
dent of position angle. Although this assumption might 
not be verified for single galaxies, it is always automat- 
ically realized in the case of stacked images. Let us 
now consider a set of QSOs homogeneously distributed 
in sight-lines around L galaxies. The number of QSOs 
that arc actually absorbed at given impact parameter R 
is thus given by 

dNabs {R,L)^ Nabs (L) dP{R, L) (6) 
^2TTNabs{L)p{R,L) RdR, (7) 

where Nabs{L) is the total number of QSOs with detected 
absorption. If we now stack all the absorbed lines-of- 
sight, we find the azimuthally averaged surface bright- 



2^ R dR-NabsW 

This is valid in the limit that the effective extent of each 
galaxy is significantly smaller than R. As galaxies with 
L ^ L* have typical effective radii of the order of 3 kpc, 
at i? > 10 kpc Equation [5] is a fair approximation. Sub- 
stituting Equation [7] into this expression we see that the 
area from the probability cancels out with the normal- 
ization of the azimuthal average and we have 

SB{R,L)^L-p{R,L). (9) 
Therefore, the SB profiles presented in the previous fig- 
ures are directly proportional to the probability of pro- 
ducing absorption at a given distance from a galaxy. If 
we now consider contributions from galaxies of different 
luminosities, we can write 

SB{R) = jL dL, (10) 

that is to say that the SB at any radius R is proportional 
to the luminosity-weighted impact parameter distribu- 
tion. 
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Fig. 8. — The correlation between the shape of the surface brightness profiles and the absorption strength. In the top panels, the slope 
of the best fit powerlaw slope (20 kpc < R < 100 kpc) is plotted against the absorption strength, for the three z^bs bins, as indicate by the 
labels in the middle. Blue open squares, green filled circles, orange empty stars, and red filled triangles are shown in the upper row for the 
g-, r-, and z-band respectively. The black filled diamond in the lower row is the weighted average of the four bands. The horizontal error 
bars show the extent of each EW bin. The bottom panels are the same as the top ones, but for the first moment of the SB distribution Ri 
(see text). 



In order to take into account the effect of satellite 
galaxies, or environment if a Mgll absorber is associated 
with more than one galaxy (for example in a group), 
we can extend the argument that leads to Equation [9l 
Let now L be the luminosity of the brightest associated 
galaxy and R measure the impact parameter from it. The 
presence of secondary galaxies increases the total lumi- 
nosity per absorber and produces a spatially extended 
contribution to the average SB. By defining Ki{R) as 
t_he average surface brightness around the location of an 
L galaxy, including all its satellites and possible compan- 
ions, and normalized dividing by L, we can then rewrite 
Equation [5] as 

SB{R,L) = Ci-L-{p{R,L)®Ki{R)) (11) 

where the Ci accounts for the increased total luminos- 
ity per absorber and p{L,R) is convolved with Ki{R). 
Finally, integrating over luminosity we obtain the analo- 



gous of Equation [To] for multiple associated galaxies 
d{p{R, L)®Kl{R)) 



SB{R) = J Cl-L- 



dL 



dL 



(12) 



It is immediately seen that this last equation simplifies 
to Equation [TUl when Cl = 1 and Kl{R) = Sj:,irac{R) for 
every L, i.e., when there is only one associated galaxy 
per absorber, with negligible spatial extent. 

4.2.1. How many galaxies per absorber? 

Although the association of an individual galaxy with 
an absorber is a somewhat ill-defined procedure, in or- 
der to interpret the SB profiles and especially the SEDs, 
it is highly relevant to understand how many (bright) 
galaxies per absorber contribute to the cross-correlating 
light. Preliminary results by Menard et al. (in prepa- 
ration) indicate that at Zabs < 0.6 there is on aver- 
age one bright galaxy detected in excess within ~ 100 
kpc around absorbed QSOs in comparison to reference 
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QSOs. As we will show in ij5l the average total lumi- 
nosity < Cl ■ L > (within 100 kpc, in the \ow-Zabs bin) 
is rrir = —21.5 (±0.2) or « 2.8 x 10^°^©,^, correspond- 
ing to th e L* of the local field luminosity function in the 
SPSS by'Blanton et al.' (2003) (we adopt iBlanton et aH 
l2005l to convert from magnitudes in the z=0.1 system to 
the rest- frame). If we compare this value to the LF at 
z ~ 0.5, then thi s luminosity corr esponds to 0.5 L* (e.g. 
in COMBO-17, IWolf et al.li200l . Such a relatively low 
luminosity inside a 100 kpc circle shows that, on aver- 
age, strong Mgll absorbers with O.SA < Wq ^ 3^ do 
not reside in dense galaxy environments. This statement 
however cannot be proved for stronger systems (with 
Wo > 3A) alone as their relative contribution to the sig- 
nal is substantially lower. 

Therefore, in the following we will interpret the SB 
profile SB{R) to first order approximation as the lumi- 
nosity weighted average probability of intercepting an 
absorber at a given impact parameter R from a bright 
galaxy. Regarding the corresponding spectral energy dis- 
tribution (SED), the inclusion of small companions and 
satellite galaxies in the integrated light is expected to 
have a minor influence on our estimates of total luminos- 
ity and spectral properties. In practice we will treat such 
properties as the properties of a single absorbing galaxy. 

4.3. The impact parameter distribution of Mgll 
absorbers 

Based on the arguments presented in the previous 
paragraphs, we translate the observed SB distributions 
around the absorbed QSOs into a luminosity weighted 
probability of intercepting an Mgll-absorbing system at 
a given impact parameter from a galaxy. In differential 
form this is directly given by the SB profiles (see 
modulo an arbitrary normalization factor. As already 
noted, simple power laws are able reproduce the SB pro- 
files, and hence also the differential probability profiles, 
reasonably well. Here we introduce the integral form 
R 

P{< R) = J 27r R' dR' J L^^^^^dL (13) 



where P{< R) is the luminosity weighted probability of 
measuring an impact parameter between and R. This 
form of cumulative probability distribution is particu- 
larly useful for defining characteristic scale lengths. Ac- 
cording to EguationslSland llOi this is simply proportional 
to the apparent luminosity integrated between and R, 
F{< R). The proportionality factor between the two 
distributions is just given by the total integral flux. 

Observationally, however, we are faced to two limita- 
tions. At large radii small background uncertainties yield 
huge uncertainties in the estimate of the total integral 
flux. Therefore we are forced to limit the flux integral to 
an arbitrary large radius where the background level un- 
certainty is negligible. This is chosen as 100 kpc, which 
implies that the impact parameter distribution will be 
normalized to 1 between and 100 kpc. The second lim- 
itation is that at small radii {R < 10-20 kpc) the finite 
size of the PSF does not allow a reliable estimate of the 
flux (see i i4.ip , and therefore we can produce flux growth 
curves only from 10-20 kpc outwards. 

We overcome this problem by extrapolating the growth 
curve to i? = using an analytic fit to the empirical 



13 

growth curve. Only the conservative range between 20 
and 100 kpc is used. As a convenient fitting function we 
choose an hyperbolic tangent generalized as follows: 

F{R) = a tanh[/3 • (i? - ,5)] -f 7 (14) 

From this fit we can derive the (extrapolated) total flux 
between and 100 kpc Fmo = F(IOO) - F(0), and hence 
express 

P(<i?) = citanh[c2-(i?-C3)]+C4, (15) 

where 

Cl = a/i^ioo (16) 
C2 = /3 (17) 
C3 = 5 (18) 
C4 = Cl • tanh(c2 • C3). (19) 

The i-band photometry in the low-redshift sample is 
used to derive the cumulative impact parameter distri- 
bution P(< R). Thus what we present here is weighted 
on the rest- frame g-band light. In table Owe report the 
four coefhcients of the analytic formula given in Equa- 
tion[T5]for several bins of Wq- These analytic expressions 
are meant to be handy representations for the P(< R) 
extracted from a robust statistical analysis. No physical 
model is implied. However, the reader should be aware 
that below 20 kpc these formulae are extrapolated and 
that the normalization within 100 kpc is arbitrary. 

In the left panel of Figure [9l P(< R) is plotted for 
all EWs as a solid line. The dark grey shaded area is 
the 1-cr confidence region for P(< i?) as it is directly de- 
rived from the i-band (roughly rest-frame (7-band) flux 
growth curve of the entire sample at low redshift. The 
integral impact parameter distributions for the low- and 
high- Wo systems are overplotted in the same panel with 
dashed red and dotted blue lines, respectively. Again, 
the higher concentration of stronger systems is evident. 
This is also quantified by the values of P50 and -R90 re- 
ported in Table [2] that represent the distances from the 
"parent" galaxy within which the probability of finding 
an absorber system is 50 and 90% respectively (as com- 
puted from Equation I15p . 

In the right panel of Figure[9l we differentiate P(< R) 
to derive the differential impact parameter distribution 
^(i?). As in the left panel, the different lines repre- 
sent the three EW samples. The dark grey shaded area 
in this plot shows the normalized SB profile (1-cr confi- 
dence level) for the "AU-EW" sample. The agreement 
with the corresponding analytic curve (black solid line) 
is remarkable and extends even beyond the fitted range; 
this testifies to the goodness of the fit. 

Three lines of different logarithmic slopes are also plot- 
ted in figure ini (right panel); they can be compared with 
the lines representing the analytic tanh fit, and with the 
power-law best fit slopes reported in Table [1] column 4. 

It is instructive to compare our results with previous 
studies on Mgll impact parameters. For thi s, the main 
sampl e available in the literature is that of iSteidel et al.l 
(| 19941 ) with about 50 "identified" Mgll absorbing galax- 
ies of whi ch 70% have spectroscopic redshifts. As de- 
scribed in ISteidel (l995), the global picture that arose 
from their study implies that bright galaxies are sur- 
rounded by a halo of gas detectable through strong Mgll 
absorption lines up to about 40 kpc. More specifically. 
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Fig. 9. — Luminosity-weighted impact parameter distributions for the low Zabs sample, normalized between and 100 kpc. The lef t panel 
shows the cumulative probability of observing a galaxy at impact parameter < R from an Mgll-absorbed QSO (see Equation I13I I. The 
three lines are the tanh analytic fits to the distributions for all (black solid line), high-Wo (blue dotted line), and low-lVo (red dashed line) 
systems. The dark grey shaded area represents the l-cr confidence interval empirically derived from the photometry for the whole low-z^j^, 
sample. Data points in the regions shaded in light grey were excluded from the fit. The right panel reports the differential probability 
^(ij), as analytically derived from the cumulative one. Different line styles and shadings represent the same quantities as in the left panel. 
Three lines with diff'erent logarithmic slopes are also plotted for easy comparison. 



TABLE 2 

Parameters of the impact parameter distribution P{< R) 
(0.37 < Zabs < 0.55) 



Wo(A2796) 


Cl 


C2 


C3 


C4 


^?50 


R90 






kpc 


-1 


kpc 




kpc 


kpc 


All 


7.900 


7.275 • 


10-3 


-163.91 


-6.567 


34.6 


81.4 


LOW-VKo 


13.555 


7.406 • 


10-4 


-29.98 


-0.301 


49.9 


90.0 


HIGH-Wo 


11.177 


1.359 ■ 


10-2 


-109.89 


-10.103 


23.9 


68.5 



using several quasa r sightli n es wi thout any strong Mgll 
absorption feature, ISteidell ()1995( ) derived that the size 
of Mgll halos follows R{Lk) ^ 38 h'^ (LK/Ll^f-^^ kpc. 
It is interesting to note that the mean impact parameters 
inferred from his study are similar (although systemati- 
cally lower) than the ones obtained in our analysis. 

However, whereas Steidel et al. did not report detec- 
tions of galaxy impact parameters larger than 40 
kpc, we find a significant SB signal up to about 200 kpc. 
More specifically we find a similar amount of correlated 
light on scales r < 50 k pc and 50 < r < 100 k pc. As re- 
cently pointed out by Churchill et al] (|2005af) the study 
performed by Steidel et al. suffered from several selec- 
tion biases. For example, once a galaxy was identified 
at the absorber redshift, no further redshifts were deter- 
mined in the field thus discarding the possible detection 
of several galaxies contributing to the absorption. More- 
over galaxies were targeted depending on their colors and 
magnitudes, and spectroscopy was performed beginning 
with the smallest impact parameters. Churchill et al. 
have started to reinvestigate some of these absorber sys- 
tems more systematically. They have reported several 
mis-identifications. Their preliminary results show the 
existence of more than one galaxy at the absorber red- 
shift in some cases, as well as a number of impact pa- 
rameters as large as 80 kpc. 



In light of these results, it is worth emphasizing the 
power of our stacking approach and its ability in pro- 
ducing unbiased results over a range of impact parame- 
ters that are hardly accessible with any other "classical" 
method. Our analysis provides us with a well-defined 
cross-correlation between gas and light: {Wo{(t)) L{(t)+9)). 
Trying to identify a galaxy responsible for the absorp- 
tion is not awell-defined process. On the one hand, there 
can very well be galaxies at smaller impact parameters 
that are too faint to be detected. On the other, Mgll 
clouds could be at such large distances from the associ- 
ated galaxy that confusion can arise due to galaxy clus- 
tering. Correlation functions (using a stacking technique 
or galaxy number counts) are needed in order to probe 
the statistical properties of connections between absorp- 
tion and emission properties. 

5. INTEGRAL PHOTOMETRIC PROPERTIES OF MGII 
ABSORBING GALAXIES 

In this section we focus on the integrated absorbers' 
light, with the aim of characterizing the absorber-related 
galaxies in terms of spectral energy distribution (SED) 
and total optical luminosity. Our approach consists of 
comparing the observed fluxes with suitably transformed 
galaxy templates observed in the local Universe. Our ap- 
proach and the subsequent interpretation are based on 
the assumption that there are galaxies in the local Uni- 
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verse whose SED closely reserables the average SED of 
the galaxies connected to the Mgll absorbers at ^; > 0.4. 
Moreover we assume that evolution within a redshift in- 
terval Ax = 0.2 is negligible (which will be verified as 
shown below). We begin by selecting 39 nearby galax- 
ies spanning the whole range of morphology, star forma- 
tion history and current activity, for which good quality 
spectral coverage from t he UV to the near IR is pro - 
vided by [Calzctti et all (Il994j). iMcQuade et al] (|1995D . 
and IStorchi-Bergmann et al.l ( 19951 )^". Let us now con- 
sider a given sample of N absorbers and its redshift dis- 
tribution {zabs,i}i=i,N- For each template SED, we com- 
pute the effective SED which would be observed if all 
absorbing galaxies had this same template spectrum and 
the same absolute magnitude in the rest-frame. 

The expected effective flux density at the observed 
wavelength A is then given by: 



/off (A) 



10 



-0.4AM 



N 



N 



i=l 



1 ~(~ ^absA 



(20) 

where /o is the flux density of the template observed 
at the canonical distance of 10 pc with a fixed absolute 
magnitude Mg, N is the number of absorbers in the sam- 
ple, fi{zabs,i) is the distance modulus corresponding the 
given absorber redshift, and AM is the difference be- 
tween the actual absolute magnitude of the absorbing 
galaxies (which is assumed to be the same for the all ab- 
sorbers in the sub-sample) and Mq. AM is treated as 
a free normalization parameter in the following. From 
the effective SED, synthetic observed- frame flux densi- 
ties in the four bands are computed and compared to 
the observed fluxes from the stacking, using the stan- 
dard definition of for correlated errors: 



X 



(21) 



b—g,r,i,z b'—g,r,i,z 



where A/b are the differences between the observed and 
the effective template fluxes, and ab,b' is the inverse co- 
variance matrix of the observed fluxes in the four bands. 
For each template, we determine the normalization factor 
AM that minimizes the x^, and consider this minimum 
as a measure of the goodness of the flt attainable with 
that given template. We illustrate the main re- 

sults of this analysis in Figure [TOl Here we consider the 
three Zabs subsamples. In the three panels in the top row 
the observed flux densities, derived from the integrated 
magnitudes between 10 and 100 kpc (see table [3]) in the 
three Zabs bins, are reported as the orange filled dots with 
error-bars (the horizontal error-bars represent the effec- 
tive width of each passband). For the sake of illustration, 
we now consider only four SED templates that are repre- 
sentive of the whole range of galaxies. Their rest-frame 
spectra are shown in the four panels of the bottom row. 
From left to the right, we consider: 

a) SL typical early- type spectrum (from the SO galaxy 
NGC 1553), characterized by a very suppressed UV flux, 
strong 4000 A break and metal lines, and flat shape at 
A > 6000 A. 

The observed spectra are shifted back to their rest frame 
and corrected for foreground Galactic e xtinction using the 
ICardeUi et al.l I I1989I ') extinction curve and the lSchlegel et al.l lil998l ') 
dust distribution. 



b ) the spectrum of an intermediate spiral galaxy (from 
the Sb galaxy NGC 7083), which is characterized by 
stronger UV flux, a weaker 4000 A break, weak emission 
lines, and very moderate blue slope in the continuum at 
A > 6000 A. 

an early type spiral (NGC7552, SBab) which is 
undergoing a starburst, partially attenuated by dus t 
(StarBurst class 6 in the atlas of 'Calze tti et al.lll994f ): 
the unattenuated blue continuum longward of the 4000 
A break and the strong emission lines coming from 
the star-burst coexist with a strong 4000 A break and 
strongly suppressed UV flux, typical of the unextincted 
old stellar population. 

d) Finally, the fourth panel shows the spectrum of an 
unextincted star burst galaxy (N GC 5253, starburst class 
1 in the atlas of ICalzetti et al.l[l994 ). which is charac- 
terized by blue optical continuum, strong emission lines, 
and highly enhanced UV flux. 

These four spectra are redshifted and convolved with 
the distribution of Zabs and normalized to the observed 
data points as explained above. The results are reported 
in the three top panels of Figure [TOl where different line 
styles and colors code the four original spectra. We also 
report the relative to each model. It is apparent that 
the two extreme templates (passive early type and pure 
star burst) completely fail to reproduce the observations 
at any redshift. Good fits are obtained only with the 
intermediate types, with the best results coming from 
the intermediate spiral with partially obscured starburst. 
This latter young component appears to be particularly 
required in order to reproduce the observed "jump" due 
to the 4000 A break, which in addition, is required to 
be weaker than in a normal Sb galaxy. A striking result 
from Figure [TO] is also the fact that the same template is 
the best fit to the observed fiuxes at all redshifts. In other 
words, we do not observe any significant evolution in the 
average SED of galaxies associated with Mgll absorbers 
from z = 0.4 to z = 1, which validates our assumption 
of negligible redshift evolution. 

From Figure [TO] we also note that the SDSS bands 
mainly constrain the rest-frame u and g flux. In the 
highest redshift bin the situation is even more extreme, 
as the g-, r-, and i-bands sample the flux blue- wards of 
the 4000 A break, while the z-band is not very sensitive; 
this results in a very poor constraint to the global proper- 
ties of the stellar populations, and in an over-sensitivity 
to the effects of dust attenuation in the rest-frame UV 
at Zabs > 0.6. 

5.1. SED fitting analysis 

We now extend our analysis, using the complete library 
of the 39 UV-optical spectra selected as reported above, 
to assess the previous result statistically, and in terms 
of rest-frame colors and luminosities. All templates are 
classified by morphology (according to NED^^) and their 
rest-frame colors, u — g, g — r, r — i^ are synthesized from 
the spectrum. The distribution of the templates in this 
3-dimensional color space is plotted in Figures [TT] and 
[TOl for different absorber subsamples, respectively. In 
each of these figures, the morphological type of the SED 
template, from EUiptical-SO to Irregular/Blue Compact 
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Fig. 10. — Illustration of the SED fitting technique. In the bottom row we plot the rest-frame spectra (UV to near-IR) of four galaxies 
which are representative of the whole range of galaxy SEDs, from the early types (to the left) to the late and star-bursting types (to the 
right). The flux normalization is arbitrary. In the top panels the filled dots with error bars display the flux density in the four SDSS bands 
(observed frame) for the three z^bs bins (low is to the left and high to the right). The horizontal error bars span the filter band- width. 
The lines represent the best-fit renormalized spectra computed by convolving the redshifted templates of the bottom row with the redshift 
distribution of the absorbers (see text for details). The line styles and colors are the same as in the bottom panels. The of each model 
is reported as well. Note how much better the intermediate type templates (NGC 7552 in particular) match the observational data at all 
redshifts. 



Dwarf galaxies, is coded according to the symbol shape 
(from circles to triangles) and color, whereas the size 
scales inversely to the x^i as indicated in the legend. The 
best-fit model is identified by the cross-hairs; all models 
that result in < 6 (or a per degree of freedom less 
than 2) are identified by the open symbols, while other 
models are represented by the filled ones. We im- 
mediately observe that the rest-frame u — g color, and, 
to a lesser extent, the g — r are well constrained by our 
observations, and indicate that the average SED of the 
absorbing galaxies is consistent with intermediate-type 
spiral galaxies (Sab to Sbc) of the local Universe at all 
studied redshifts (Figure fTTjl. We consider the colors of 
the best-fit model as best estimates of the rest-frame col- 
ors and = 6 as a confidence limit, and report these 
values in columns 2, 3, and 4 of Table Uf^. We conclude 

More specifically, the highest value corresponds to the reddest 
rest-frame color of any model for which < 6. Vice versa, the 
lowest value is for the bluest color of any model with < 6. Note 
that the highest (lowest) values for two different colors do not need 



that in the entire redshift range 0.36< z <1.0 the av- 
erage SED of the absorbing galaxies is characterized by 
rest-frame colors u — g — 1.1 ± 0.2, g ~r — 0.8 ± 0.2, and 
r -i ^0.3 ±0.1. 

By applying the same kind of analysis to the two sub- 
samples of low-tVo (Figure [12] left panel) and high-Wo 
systems (Figure fT2l right panel) at low redshift, we ob- 
serve significant differences between the two classes; i.e., 
the colors of low-M/p systems are best fitted by the tem- 
plate of an elliptical galaxy and the high-Wo systems 
appear to be significantly bluer. This is shown quanti- 
tatively in Table 2] where we can see strong differences 
in u — g between the two subsamples. At low redshift 
high-Wo systems are about 0.6 mag bluer than low-Wo 
ones. At higher redshift the signal is too weak to perform 
a similar analysis. 

The difference between the SEDs of low- and high- 
Wq systems is particularly well illustrated in Figure [T51 

to be taken from the same template. 
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TABLE 3 
Integral photometry and colors 



VKo(2796) bin 


Radial range 


g 


Integral flux (mag) 
r i 


z 


g_r 


Color (mag) 
r-i 


i-z 










0.37 < 


: Zabs < 0.55 










Low-VKo 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


23 50+°1^ 
23.88_Q 


22 00'^?:^,i 

''^•""—0.17 
91 01 -1-0.30 

2i.i4«:- 


21 40-^°-?J. 
21 12+0-29 

20.50l0;22 


20 82-^0-29 

90 59+0.53 
ZU.OZ_o 3g 

19 01+0. 36 


1.50 ± 0.34 
2.08 ± 0.99 
1.78 ±0.53 


0.60 ± 0.12 
0.69 ±0.18 
0.65 ±0.13 


0.57 ± 0.23 
0.60 ±0.39 
0.59 ±0.29 


Intermediate- VFo 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


22.95+^-j» 
24.38l-o'4 


22.18--1; 

99 Vl ~t"0.73 
' -'--0.43 

21.66l°;i 


21.88t"oi^ 

99 97+0.96 
ZZ.Z(_Q 50 

21 30+0-41 

Z1.0U_Q 3Q 


22 02-'-^ '^'' 
ZZ.UZ_o 57 

91 95+1-15 
Z1.Z0_Q 55 

20.811°;!^ 


0.76 ± 0.22 
1.67 ± 1.60 
1.03 ± 0.44 


0.31 ±0.19 
0.44 ± 0.48 
0.36 ±0.26 


-0.14 ±0.67 
1.02 ±0.78 
0.49 ±0.61 


High- Wo 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


22 54+"-^'-^ 

99 1 9+0.35 


21.58_o--^3 

99 1 o-t-0.41 

^^.10_Q 3Q 


21 09-'-"-^^ 
zi.ua_Q j^3 

91 85+0.54 
Z1.00_o 3g 

20.65-'-oi7 


20.59^^:^8 
21 35+1-24 

Z1.00_o 5g 

20 15+0-39 
ZU.iO_o 29 


0.97 ±0.14 
1.21 ± 0.59 
1.05 ± 0.25 


0.48 ± 0.08 
0.28 ±0.31 
0.41 ±0.14 


0.50 ±0.18 
0.50 ±0.70 
0.50 ±0.31 


All 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


22.89+0;i^ 

oo '7c;-t-0.82 

99 4Q+0-27 
zz.'±y_Q 21 


22.16l°:g 
21 25-'-°-12 


21 39+0-11 
21-67l°:^5 
20.7710:1^ 


20.97t^;i« 

21 m+°-*° 

Zi.UU_Q 29 

20 93+0-24 
zu.zo_o 20 


1.03 ± 0.11 
1.59 ± 0.50 
1.24 ±0.21 


0.47 ±0.07 
0.49 ±0.16 
0.48 ±0.09 


0.42 ±0.15 
0.67 ±0.32 
0.54 ±0.21 


0.55 < Zabs < 0.76 


low-VFq 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


24 07+1-31 
^^•'^'-0.58 

23.561°;^? 


23.48l0:« 
23.6llJ:^^ 

99 7Q+0.59 

3g 


99 Q-3+0.41 

zz. oo_0 3Q 
22.43l°;| 
21.86lO;^8 


22 fiO+l l* 

zz. DU_Q 55 

91 1Q+0.52 
zi.ia_o 35 

20 93+0-54 
zu.»o_Q 3g 


1.16 ± 0.60 
0.45 ± 0.80 
0.77 ±0.59 


0.65 ±0.29 
1.18 ±0.59 
0.93 ±0.35 


0.23 ±0.70 
1.24 ±0.51 
0.93 ±0.47 


Intermediate- Wo 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


94 14+O.M 

24 •^n+2-40 
23.46i°;^^ 


22.90+Qig 
23.73lj;^6 

22.49tQ;2? 


99 9O"l-0.23 
ZZ.ZO_o xg 

99 +0.56 
zz. 01_o 37 

2i.64i°;g 


22.02+o-i^i 
22.34l^--20 

21-421^6 


1.23 ± 0.38 
0.56 ± 1.01 
0.97 ±0.57 


0.62 ±0.17 
1.22 ±0.59 
0.85 ±0.25 


0.26 ±0.39 
0.17 ± 1.09 
0.22 ±0.60 


High- Wo 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


OO otr + U.^i 

23.7810-^^ 
22.79l°;^i 


22.33loii 

99 71 +0.35 
' -^-0.26 

91 75+0.17 
'O-0.15 


21.81±li:}^ 

22.08-'-o27 
21 19+0-18 

Zl.iM p jg 


21 74+0-^'/ 

22 39+ '"1 
zz.oy_o 79 

91 97+0.73 
-0.43 


1.02 ± 0.17 
1.08 ± 0.47 
1.05 ± 0.25 


0.51 ±0.10 
0.63 ±0.26 
0.56 ±0.14 


0.07 ±0.30 
-0.31 ± 1.11 
-0.08 ±0.52 


All 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


23.S»_0 

94 m+0-50 

95 1Q+0-29 


22.79+'^;ii 
9Q 91 -t-0.32 
zo.zi_o 25 

99 9Q + 0.16 
ZZ.ZO_Q j4 


99 99H-0. ^2 
ZZ.ZZ_o j^Q 

99 ,1+0.24 

zz. 01_Q jg 

21 51+0-14 

Z1.0±_Q J3 


22.06+0 22 
91 HI +0.50 

91 1 '7+0.33 


1.10 ±0.17 
0.80 ±0.39 
0.97 ±0.23 


0.57 ± 0.09 
0.90 ±0.22 
0.72 ±0.12 


0.16 ±0.24 
0.50 ±0.41 
0.34 ±0.29 








0.76 < 


: Zabs < 1-00 










LOW-Wq 


10- 50 kpc 
50-100 kpc 
10-100 kpc 


24 2fi+°-'*^ 
24.87l^"4 
23.771°;^° 


24 53+1-10 
24 95+ '"1 

90 Qr'-t-2.30 

ZO.»l_Q gg 


90 90+0.41 
ZO.ZO_Q 3(j 

23.03l°;P 


22.65«:S 

23 1 1 + '"1 
zo. 02 

22 lO-'-l '''^ 

zz. iU_og4 


-0.27 ±0.66 
-0.08 ± 1.80 
-0.20 ±0.95 
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10- 50 kpc 
50-100 kpc 
10-100 kpc 
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0.24 ±0.21 
0.57 ±0.71 
0.36 ± 0.33 


0.82 ±0.18 
0.94 ± 0.48 
0.87 ±0.25 


0.53 ±0.29 
0.84 ±0.65 
0.68 ± 0.40 



where we plot the best fitting SEDs for the three bins 
of Wo at low redshift. The top half of each panel repro- 
duces the observed photometric data points overlaid on 
the convolved and normalized template spectrum, while 
the bottom half displays the original template spectrum. 
It is apparent that the weakest systems are best repro- 
duced by a red passive template, while at Wq > l.lA the 
SED closely resembles those of local actively star-forming 
galaxies. 

From the fitting of the SED we also derive the normal- 
ization factor, which enables us to compute rest-frame 
luminosities in different bands. As mentioned above, the 
best constraints are obtained for the rest-frame u— and 
g— band fiuxes, which correspond to the central and most 



sensitive r- and i-band in the observed frame. The uncer- 
tainty in these fluxes is driven mainly by the photomet- 
ric errors. On the other hand, the rest-frame r— and i— 
band fluxes are extrapolated with the aid of the template 
SED, and therefore are subject to the additional uncer- 
tainty of the choice of the SED template. This is clearly 
the major source of error in the reddest rest-frame bands 
at high redshift. In Table IH we report the rest-frame ab- 
solute magnitudes in u, g, r, and i in columns 5, 6, 7, and 
8, respectively, for different absorber subsamples. The 
uncertainty in the normalization factor, which is given 
purely by the photometric error, is given in column 9. 
This is computed as the variation of the normalization 
factor which is required to increase the value of by 
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Fig. 11. — "Goodness of fit" distribution in the rest-frame color space for the three z^^a 
galaxy spectrum (taken from fCalzetti et al.lll993 : IMcOuade et al.lll995l : IStorchi-Bergmann et alJIiggSI I. whose corresponding morphology 
is coded by its shape and color, as indicated in the legend. The size of the symbol is inversely related to the of the template, as from 
legend. The best fitting template is marked by the cross-hairs, while all templates yielding < 6 (our confidence limit, see in the text for 
details) are coded with open symbols. All other templates arc shown as filled symbols. 



1 with respect to the minimum while adopting the best 
fit template. All the following considerations are made 
using photometric quantities integrated between 10 and 
100 kpc. Looking at the rest-frame luminosity in u and g 
as a function of rcdshift we note a significant trend with 
higher redshift systems having higher luminosity; the in- 
crease is roughly 30 to 50% from z ^ 0.4 to z ~ 0.9. 
It must be noted that such a luminosity evolution has a 
negligible impact on our interpretation of colors within 
each Zabs bin; in fact, taking luminosity evolution into 
account changes the observed-frame colors synthesized 
from the model SEDs by 0.02 mag at most for interme- 
diate galaxy spectral types. It is worth noting here that 
luminosities and colors are scale-dependent. This effect 
is illustrated in Table [4] where we report the colors and 
luminosities for the low-Zaf,., sample measured between 
10 and 50 kpc. With respect to the total quantities the 
rest-frame luminosities are lowered by 0.45 (u-band) to 



0.83 mag (i-band), while the colors are systematically 
bluer. This color gradient is primarily a consequence of 
the fact that strong systems, which are also bluer, are 
more concentrated and dominate the total light in the 
inner tens of kpc. 

5.2. The SED and luminosity of Mgll absorbing 
galaxies in context 

At the beginning of this section we presented an in- 
terpretation of the colors of the absorbing galaxies in 
terms of local galaxy spectral templates. We showed that 
the average SED is typical of intermediate type galax- 
ies, with both old stellar populations and a starburst 
component. This mix varies going from low- to high-Wo 
absorbers, with the latter having an enhanced starburst 
compon ent. It is interes ting to note that, in a different 
context. iKacprzak et al.i (2006) have recently obtained a 
result that goes in this same direction, showing a posi- 
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TABLE 4 
Rest-frame colors and luminosities 



Wo (2796) bin u — g g — f r — i u g r i 5M 

mag mag mag mag mag mag mag mag 
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five correlation between Wo and the ratio between the 
asymmetry parameter, A, and impact parameter for a 
sample of 24 intermediate redshift absorbing galaxies. 
By interpreting A as an index of galaxy "activity" , the 
observed correlation is fully consistent with our finding 
that higher Wq absorbers are associated with galaxies 
having a smaller impact parameter and a more "active" 
SED (in the sense of star formation). The correlation 
between star formation and E W is further s uppor ted by 
the observation, reported by iNestor et al.l (|2005f) . that 
strong systems occur with a relatively higher frequency 
at high redshift > 1, thus paralleling the evolution of the 
cosmic star formation rate. 

No evidence is found for significant redshift evolution 
of the average SED. This is an argument in favor of the 
thesis that Mgll absorption systems trace a well defined 
evolutionary phase of galaxies, rather than just a partic- 
ular class of halos or galaxies in a fixed mass range. 

In order to characterize this evolutionary phase in 
the cosmological context, we compare the average col- 



ors and luminosities of the absorbing galaxies with a 
well-understood sample of / + i?-band selected galax- 
ies from th e FORS Deep Field survey (hereafter FDF, 
iHeidt et~aT., 2003). In addition to multi-band photometry 
from the U to the ifg-band, the FDF catalogs provide the 
photometric redshi ft of each gala.xy an d the best fitting 
SED template (see lBender et alll2001|) . These can then 
be used to derive the rest-frame photometric properties 
of the galaxy, and to convert the observed magnitudes 
into the SDSS photometric system. We select all galax- 
ies with absolute rest-frame i-band magnitudes brighter 
than —17.0, and split the sample into the same three red- 
shift bins adopted for the absorbers. Figure [HI shows 
the oh served- frame color-magnitude diagram (CMD) of 
the FDF galaxies in the three redshift bins. Observed- 
frame bands and colors for different redshifts are chosen 
so as to correspond to similar rest-frame spectral regions, 
namely, the magnitude on the x-axis samples the flux at 
~ 45OOA, while the color on the y-axis roughly corre- 
sponds to rest-frame u — g in all three panels. The size of 
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Fig. 13. — The best fitting templates for the three Wo bins at low redshift, from the weakest systems (to the left) to the strongest ones 
(to the right). In each graph, the rest-frame template spectrum is plotted in the bottom panel, with some arbitrary intensity normalization, 
and identified by the NGC name of the galaxy. The top panels show the measured photometric points (red diamonds, whose vertical extent 
indicates the absolute photometric error in each band individually) overlaid on the convolved and re-normalized template spectrum (see 
text for details). The blue points with error bars show the error relative to the r-band data points. The shaded area is the l-a confidence 
range allowed for the normalization. The best fitting rest- frame g magnitude is also reported, along with the (non-reduced) x^- 



the square is indicative of the absolute g-band luminosity 
of each FDF galaxy (as shown in the figure legend), while 
the color of the symbol codes the rest-frame u — r, which 
is particularl y suited to separa te late- from early-type 
galaxies fe.g. lBaldrv et al.ll2004[ ). The photometric mea- 
surements relative to our absorbers, including all Wq, are 
indicated by the grey shaded ellipses in each panel, with 
the shaded area representing the l-cr confidence limits. 
For the low- z^bs bin (upper left panel) we also ovcrplot 
the two ellipses relative to the high-Wo systems (solid 
line) and the Iow-Wq ones (dashed line). Note that we 
always refer to photometric quantities within the 100 kpc 
aperture, unless specified otherwise. 

It appears that the average colors of the absorbing 
galaxies are intermediate between the very red ones 
{{u — r)rf > 1.6) of passive galaxies on the red sequence 
(the ridge of yellow-red points visible in the two top pan- 
els of Figure [H]) and the bluer colors, that characterize 
star forming spirals of intermediate luminosity. This, 
again, is seen at all redshifts. We can thus conclude that 
the average colors of the absorbing galaxies are interme- 
diate not only in absolute sense (as we already saw from 
the fit of local SED templates), but also relative to the 
galaxies of similar luminosity at the same redshift. Fig- 
ure [T3] also supports the idea that Mgll absorbers reside 
in halos hosting a, L < L* galaxy. In fact, if the total 
observed luminosity arose from a large number of low- 
luminosity galaxies, then we would expect much bluer 
colors than observed. Conversely, we note that the es- 
timated absolute magnitude of the absorbing galaxies 
{Mi^rf ~ —21.7) is perfectly consistent with those of 
galaxies sharing the same region of the color-magnitude 
diagram^^. 

Figure 1141 and the low-Zabs panel in particular, at 
first glance suggests that the absorbing galaxies could 
be those rare galaxies that populate the so-called "green 
valley" between the red sequence of passive galaxies and 
the blue sequence of the star-forming ones. In turn, this 

For instance, this would have not been the case if the observed 
<; — r of the low-Zaba sample was 0.7 instead of 1.25; in that case 
we would have been forced to conclude that the observed total 
luminosity is produced by several blue low-luminosity galaxies. 



would imply that the galaxies linked to Mgll absorbers 
are a kind of transition type that live in this phase for 
a very limited amount of time (on the order of a few 
hundred Myr). Alternatively, the intermediate color of 
the total sample may result from a kind of "morphologi- 
cal mix" of passive and actively star-forming galaxies. At 
the beginning of this S}5]we noted that high- Wo absorbers 
are linked to bluer and more star-forming galaxies than 
the Iow-Wq systems. In Figure [14] the location of two 
ellipses for the high- and Iow-Wq systems actually show 
that a systematic change in galaxy properties occurs as a 
function of Wq : while high- Wq systems occupy the locus 
of luminous star-forming spirals, the Iow-Wq ones overlap 
with passive red-sequence objects. In order to investigate 
whether the two classes of absorbers reflect the color bi- 
modality of the overall galaxy population we can use the 
three projections of the {g — r)-{r — — z) color space. 
Figure [T51 shows these distributions for the FDF galaxies 
in the low Zabs bin, with the same symbols and conven- 
tions used in Figure [14] and reported in the legend. As 
in Figure [HI the shaded ellipse represents the confidence 
region for the complete low-Zabs sample, while the solid 
and dashed ellipses are the confidence contours for the 
high- and Iow-Wq subsamples, respectively. In addition, 
we also report a red dotted ellipse to represent the colors 
of the total low-Zafc^ sample, but within the smaller aper- 
ture 10 < i? < 50 kpc. The {r—i)-{g—r) projection of the 
color space appears particularly suited to divide between 
red, passive galaxies and blue, active ones. The average 
colors of the total sample fall exactly in between the two 
clouds, but now stronger systems are clearly separated 
from the weaker ones: the stronger systems have col- 
ors typical of star-forming spirals, while the weaker ones 
overlap only with passive galaxies. It should be noted 
that intermediate-WQ absorbers have the same colors as 
stronger ones, thus indicating that the transition between 
active and passive types must be quite sharp and occurs 
at Wq ~ lA. This result indicates that absorbing galaxies 
are likely to be representative of a large range of normal 
galaxies with L < L* , whose color bimodality is also re- 
produced. Also, absorbers weaker than lA are most 
likely associated with galaxies that have already moved 
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Fig. 14. — The observed frame color-magnitude diagram (CMD) of galaxies selected from the FORS Deep Field (FDF), compared to our 
absorbers stack photometry in the three redshift bins. Points are FDF galaxies coded to indicate their rest-frame absolute g-band magnitude 
(size) and rest-frame u — r color (color), as shown in the legend. The shaded ellipse is the l-cr confidence interval for color /magnitude of 
our stack photometry for the "AH-" EW sample. In the low-^ats panel (upper left) we also display the fiducial estimates for the low- and 
high-Wo absorber sub-samples with the dashed and solid ellipses, respectively. 

one just considers the luminosity-weighted average color 
of absorbing galaxies within 50 kpc, instead of 100 kpc, 
bluer colors are obtained (as shown by the red ellipses in 
Figure 1151 leading to the wrong conclusion that passive 
galaxies are much less abundant among the absorbers^"*. 
Our extensive photometry represents a clear advantage in 
this sense, with respect to most of the previous spectro- 
photometric surveys that could cover only a limited range 
of impact parameters. 

Finally, we have also observed a trend of increasing 
rest-frame luminosity with redshift (30 to 50% more lu- 
minous from z ~ 0.4 to 0.9). If we interpret the presence 
of Mgll absorbers as indicative of a given evolutionary 
phase for a galaxy, this luminosity evolution can be seen 
as another phenomenon of downsizing, i.e., the occur- 



to a phase of passive evolution, while stronger systems 
mainly reside in halos where the primary galaxy is still 
blue and actively star-forming. It is interesting to note 
that the link between color bimodality and EW holds 
also when so-called weak systems (W^o < O.SA) are con- 
sidered. Using 7 absorbers with WFPC2/HST imagin; 
and spectroscopic spatial coverage iChurchill et all (120(1 
have recently shown that passive galaxies are also associ- 
ated with such weak a bsorbers. Moreover, the clustering 
analysis conducted bv lBouche et al.l (|2006f l shows a sim- 
ilar tendency, i.e., weaker systems appear to be more 
clustered than stronger ones, which suggests that they 
are related to redder galaxies. However, while the trend 
is reproduced, the large dark matter halo masses inferred 
by their analysis for the low-EW systems differ from the 
< L* luminosities found in the present study. 

It should be noted that the average colors, and hence 
their interpretation in terms of spectral types, is depen- 
dent on the choice of the impact parameter limits. If 



^■^ The origin of this bias resides in the different impact param- 
eter distributions for low- and high-iyo systems, which results in 
the latter being the dominant source at small impact parameters. 
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Fig. 15. — The three projections of FDF galaxies in the observed frame {g — r)-{ r — — z) observed color space, compared to our 
absorbers stack photometry in the low-z^j^s bin (0.37 < z^bs < 0.55). As in Figure [141 points are FDF galaxies coded to indicate their 
rest-frame absolute g-band magnitude (size) and rest-frame u — r color (color), as shown in the legend. The shaded ellipse is the l-cr 
confidence interval for colors from our stack photometry for the "All" EW sample. High- and Iow-Wq absorber sub-samples are shown with 
the solid and dashed ellipses, respectively. The red dotted ellipse is used to display the colors in the limited range 10-50 kpc. 

rence of a given evolutionary phase is delayed in smaller 
galaxies. 

6. INTERPRETATION 



In this section we attempt to interpret the observa- 
tional results reported in this paper. We briefly discuss 
possible scenarios that may reproduce the spatial and 
photometric measurements presented above, encourag- 
ing detailed modeling in these directions. 

In il4.2l we presented measurements of a spatial cross- 
correlation between Mgll absorption and light, i.e., 
galaxies, and we have shown that such a quantity is 
related to a light weighted distribution of impact pa- 
rameters. As explained above such observational results 
provide us with robust statistical constraints as they do 
not make use of any assumption about the nature, im- 
pact parameter and number of galaxies responsible for 
the absorption. Moreover, the large number of absorp- 
tion systems analyzed (about 600 to 1250 in each red- 
shift bin) allows us to measure the gas-light correlation 



up to about 200 kpc. It is important to note that the 
effects of galaxy clustering becomes increasingly impor- 
tant as we consider larger scales. The gas detected in 
absorption can be "related" to a given galaxy which can 
be closely surroun ded by a nu mber of neighbors. As re- 
cently shown by iChurchill et al.l (|2005b( ) this has led to 
several misidentifications in the past. 

Clustering effects are inherent to any statistical spatial 
studies and must be taken into account for interpreting 
results. Over the last few years, galaxy-galaxy lensing 
studies have illustrated this fact. This weak lensing tech- 
nique provides us with a cross-correlation between mass 
and light. Studies have shown that the galaxy masses in- 
ferred from this technique have a sign ificant contri bution 
from their environment ( Mandelbaurr Tet al.ll2006[ ). Sim- 
ilar effects are expected for gas-light cross-correlations. 

As we showed in i l4.21 the slope of the spatial cross- 
correlation between gas and light strongly depends on 
the strength of the Mgll absorbers. This change of slope 
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can be due either to a clustering of the gas around galax- 
ies being a function of Wo or to a change in the galactic 
environment. Detailed modeling using, for example, the 
halo-model approach will be required to accurately in- 
terpret these results. 

Regarding the photometric properties of the light 
found in the vicinity of Mgll absorbers, we have shown 
that stronger absorbers are related to bluer galaxies. We 
have found that systems with Wq > lA have colors 
matching those of star- forming galaxies with L < L*. 
Similar results based on different techniqu es were re- 
ported in previous studie s (e.g. iGu illcmin fc Bergeronl 
Il997t IChurchill et"al]ll999l ) . iBond et al.. (,20011) have also 
argued that the strong Mgll systems can arise from the 
"blow out" of gaseous super-bubbles that are produced 
by supernovae explosions in starbursting regions. Based 
on kin ematic studies of the Mgll lines, iProchter et al.l 
(|2006D have recently brought new evidence that "the 
strong Mgll phenomenon primarily arises from feed- 
back processes in relatively low mass galactic halos re- 
lated to star formation." Our results on systems with 
Wo > lA indeed provide a statistically robust support 
to their conclusion, as far as the mass and the star- 
formation activity are concerned. We have shown that 
systems with Wq > 1.5 A are found to be in the direct 
vicinity of a galaxy: 50% of the associated light is within 
25 kpc. Given typical super-wind exp ansion velocities 
of a few 100 km sec~^ ~ 1 kpc Myr~^ ()Heckman et al.l 
|2000[ ). the Mgll clouds have enough time to reach their 
observed location while the starburst is still active (as- 
suming that the duration of the burst is of the order of 
10 Myr). On the other hand we have shown that Mgll 
systems weaker than ~ 1.1 A are found at larger impact 
parameter (50% of the galaxy light is beyond 50 kpc from 
the Mgll system) , and are predominantly associated with 
reddish and passive galaxies. A possible explanation may 
be that these weaker systems are also older, i.e., observed 
a longer period of time after their outflow from the par- 
ent galaxy. These Mgll systems thus have had time to 
travel farther away from the galaxy, and in the mean- 
while the burst has finished and the galaxy has turned 
into a passive red one. The interaction of the clouds with 
the intergalactic medium might explain the lower Wq ; on 
long time-scales the Mgll clouds can be eroded by hot- 
ter intergalactic gas and the multiple components that 
give rise to large Wq are separated and diluted. Depend- 
ing on the actual duration of the burst, we can estimate 
the time-scale to turn the SED of the high-Wo absorb- 
ing galaxies into those of the weaker ones as something 
ranging from a few hundred Myr to 1 or 2 Gyr. Hence, 
the process of dilution of the Mgll clouds has to become 
effective on the same time-scale, and the expansion speed 
at distances larger than 25 kpc must progressively slow 
down to roughly 30 kpc/1 Gyr ~ 3 km sec""'^, in order 
to explain the broader distribution of impact parameters 
for weaker systems. 

In c ontrast to outflow models, iMo fc Miralda-Escudd 
(|1996l ) introduced a model in which metal absorption 
lines are related to gas infalling onto a galaxy which 
then triggers star formation. In such a scenario, weak 
systems are associated with galaxies that still have to 
accrete gas and trigger a (new) episode of star forma- 
tion. In contrast, strong systems can be associated with 



the gas that is being accreted onto a galaxy and trans- 
formed into stars. The infall model therefore explains 
the anti-correlation between EW and impact parameter 
and also the correlation between blue colors and EW. 
In this model, weak absorbing galaxies must be inter- 
preted as (the end of) a quiescent phase that precedes a 
new substantial burst of star formation, that will proba- 
bly turn the galaxy back into a star-forming spiral. The 
minimum duration of this quiescent phase can be esti- 
mated to be 1-2 Gyr, i.e., the time to turn a galaxy red 
after a substantial episode of star formation, consistent 
with the observed mean SED of weak systems. How- 
ever, such a kind of star formation history (SFH), with 
intense bursts interleaved with long quiescent periods, 
appears to be quite rare for L < L* galaxies (correspond- 
ing to the integrated light found within 100 kpc of Mgll 
absorbers). Such galaxies generally exp erience smooth 
continuous SFHs (|Kauffmann et ani2003[ ). If such an in- 
terpretation can be confirmed quantitatively, it will show 
that the simple infall model does not contain the neces- 
sary mechanisms to reproduce the observed trends be- 
tween the luminosity- weighted quantities and Wq. 

6.1. Implications for DLA systems 

Damped Lyman-a systems (DLAs) are believed to host 
most of the hydrogen in the Universe and p rovide the 
reserv oir for star formation. As shown by I Rao et all 
(|2006l ). strong Mgll absorbers are often associated with 
DLAs and the correlation between the two increases as 
a function of Wo (Mq//). It is interesting to note that 
I Rao et all ((2006') found that the mean Nhi is roughly 
constant for Mgll absorbers with Wq > 0.6 A but dra- 
matically decreases for weaker systems. This result is 
similar to the color change found between low- and high- 
Wq systems: redder galaxies are expected to have con- 
verted a significant fraction of the available neutral hy- 
drogen into stars. Moreover, while the colors of Mgll 
absorbing galaxies appear to be similar in the range 
Wq > lA, we detected a transition scale at Wq ~ lA, 
below which the colors become significantly redder. It 
would be of great interest to investigate whether other 
properties of metal absorbers or Lyman-a systems are 
seen to change around the scale of Wo{MgII) ~ 0.6 — 1 A. 
Understanding the origin of this value would give us im- 
portant insights into the nature and properties of the gas 
in and around galaxies. 

We can also note that, as the fraction of DLAs in Mgll 
selected systems increases as a function of Wo{MgII), 
our results suggest that stronger DLAs are expected to 
be associated with bluer galaxies, and more specifically, 
with more star forming galaxies. 

7. SUMMARY AND CONCLUSIONS 

In this paper we present a statistical analysis that al- 
lows us to constrain the mean photometric properties 
of Mgll absorbing galaxies and the impact parameter 
distribution of the gas. By stacking a large number of 
PSF-subtr acted images of quasars with absorbers and 
comparing them with similar but non-absorbed quasars 
we isolate the excess light associated with Mgll absorber 
systems and investigate its properties. Contrary to a 
number of previous studies which attempted to iden- 
tify individual absorbing galaxies, our method measures 
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a more general (and well defined) cross-correlation be- 
tween gas and light, and does not use any assumption 
about the extent of the gas and/or a redshift interval 
within which the galaxy can be found. Moreover, it is 
directly applicable to large samples, and can probe large 
scales. In the present study we use a sample of about 
2,800 quasar sightlines with one strong Mgll absorber 
(VKo(A2796) > 0.8 A) and a control sample of more than 
11,000 reference quasars extracted from SDSS DR4. The 
large number of systems allows us to study the proper- 
ties of the absorbing galaxies as a function of redshift and 
rest equivalent width (Wq), thus partly compensating the 
loss of information inherent to any stacking procedure. 
Our results are as follows: 

1) We obtain significant detections of cross- 
correlations between Mgll gas and light out 
to ^ 200 kpc scales. We show that the surface 
brightness (SB) distribution of such light is 
directly proportional to the impact parameter 
distribution of Mgll systems around absorbing 
galaxies, weighted by galaxy luminosity. 

2) Over the range 20-100 kpc, the SB profiles or 
impact parameter distributions are well described 
by power laws with an index a that strongly de- 
pends on the rest equivalent width of the absorbers; 
a ~ -1 for low-EW systems (O.SA < Wq < 1.12A), 
while a ~ -2 for high-EW systems {Wq > 1.58A). 
More accurate fitting formulae are also provided in 

S31 

At low redshift, i.e. 0.37 < Zabs < 0.55, where the S/N 
is high enough, we find that 

3) the mean luminosity enclosed within 10 kpc < r < 
100 kpc is Mg = -20.65 ± 0.11 for the global sam- 
ple, which is ~ 0.5 L* of the field luminosity func- 
tion in the same redshift regime. 

4) The spectral energy distribution of stronger systems 
appears to be significantly bluer. Mgll systems with 
Wo > lA are, on average, characterized by an SED 
typical of star-forming/bursting galaxies, while the 
light of weaker systems is dominated by passive 
red galaxies (although they might not dominate the 
number counts). 

Moreover, considering the entire redshift range 0.37 < 
z < 1.0, we find: 

5) no detectable evolution with redshift of both spa- 
tial distribution and global average SED, which 
is always well fitted by intermediate type spi- 
rals. This suggests that Mgll systems trace similar 
galaxies or a similar evolutionary phase of galaxies. 

6) The rest-frame optical luminosity increases by ~ 
50% from z ~ 0.4 to z < 1. This phenomenon is 
reminiscent of the so-called downsizing which is ob- 
served in a number of indicators of galaxy activity. 



We have discussed possible interpretations of these ob- 
servations including inflowing/outflowing gas and clus- 
tering effects. While detailed modeling will be required 
to robustly test each model, we point out that a sce- 
nario in which metal-enriched gas outflows from star- 
forming/bursting galaxies is in qualitative agreement 
with several types of trends seen in the data. 

The observational results presented in this paper, 
namely the cross-correlation between Mgll gas and light 
and the mean colors of absorbing galaxies, provide ro- 
bust statistical constraints of great interest for further 
modeling the gas distribution around galaxies. They will 
help produce a comprehensive and consistent picture of 
galaxy formation in both absorption and emission. 

As a side-product of this study, we have shown that 
our stacking technique is able to detect the light of QSO 
hosts and their environment, thus opening new exciting 
perspectives for the study of the relationships between 
QSOs and galaxy formation. 
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TABLE 1 

Shape parameters of the absorbers' SB profiles 
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APPENDIX 

ON THE ORIGIN OF THE SIGNAL IN REFERENCE QSOS 

In t J4.1l we noted that the PSF subtraction leaves significant SB residues even in the sample of reference QSOs. In 
this appendix we further investigate the origin of this effect and demonstrate that this is not due to a failure of the 
PSF subtraction algorithm, rather to a genuine light excess found around QSOs with respect to the actual PSF. To do 
this, we stack images of reference QSOs only, without applying any geometrical transformation to the images (apart 
from centering). We consider as possible factors of mismatch with respect to the PSF: 

(i) the redshift of the QSO. This is the case when the light excess is physically associated with the QSO, e.g. its 
host galaxy. 

(ii) the magnitude of the QSO. This is the case when cither the excess is linked to intrinsic QSO properties, or when 
it arises from a failure in deblending faint sources. 

(Hi) the color mismatch between the QSO and the star us ed to compute the PSF. In fact, redder point sources are 
expected to have a narrower PSF (see e.g. iRiedlflQBl . 

We split the three-dimensional parameter space given by zqso, the QSO apparent magnitude in i-band iqsoi and 
the color mismatch lS.{g — r) = (g — r)star — (g — i')qso into 3'^ cells including roughly an equal number of reference 
QSOs 400). Each of these 3^ subsamples is then processed and stacked separately. We take the "central" cell in the 
parameter space as a reference for the following comparisons. This is defined by 1.2 < zqso < 1-65, 18.59 < iqso < 
18.97, and 0.22 mag < A{g — r) < 0.38 mag. The central portion (30 x 30 arcsec^) of the i-band stack image for this 
reference subsample is shown in the central column of fig. [16] and repeated in the three rows. The residuals of the PSF 
subtraction are visible in these images and extend out to ~ 3 arcsec with significant SB of ~ 28.5 magarcsec"^. Let us 



z„„<1.2 z„„ai.65 
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Fig. 16. — The j-band stack images of reference QSOs in different disjoint bins in the ZQso-iQSO-^{9 ~ f") parameter space. Only the 
central 30 X 30 arcsec'^ are reproduced, using a linear grey scale to map the intensity. The corresponding levels in mag arcsec"^ are reported 
in the key bar on the right (oo means background level). The three images in the central column represent the same reference central bin 
in the parameter space, namely 1.2 < zqso < 1-65, 18.59 < iqso < 18.97, 0.22 mag < A{g — r) < 0.38 mag. The left and right panels in 
each rows are bins where only one of the three parameters at time is varied with respect to the central reference bin, as indicated by the 
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now consider subsamples where we vary the range of each of the three parameters separately, to investigate their effect 
on the residuals of PSF subtraction. In the first upper row of fig. [16] we show the images of the low- {zqso < 1-2, to 
the left) and the high- {zqso < 1-2, to the right) QSO redshift subsamples, with the same igso and A{g — r) cuts 
as in the reference. Similarly, in the second row we make the iqso vary (bright iQso < 18.59 to the left and faint 
*QSO 18.97 to the right), while keeping fix the other two parameters as in the reference subsample; and finally in the 
third row the /S.{g — r) is varied (well matched star-QSO pairs to the left, A((/ — r) < 0.22 mag, and poorly matched 
pairs to the right, A((; — r) > 0.38 mag). The immediate conclusion that can be drawn from the analysis of figHH 
is that most of the SB excess around unabsorbed QSOs is produced by low redshift QSOs. Also, high redshift QSOs 
have no SB excess at all. The most logical explanation for this is that at low redshift we do detect the light of the 
QSO's host galaxy, which drops below our detection limits at zqso ^ 1-5. The lack of any significant dependence of 
the SB excess on either A{g — r) or iqso shows that the color dependence of the PSF is well compensated by our 
algorithm (see §3.2p and that there are no serious deblending problems with fainter QSO being harder to deblend from 
neighboring faint stars or galaxies. 



